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A.  IHTRODUOTIOH 


Winds  blowing  over  the  ocean  generate  wares  on  the  water  surfaoa  whioh  are  as  vari- 
able  at  the  winds  whioh  oreate  them*  They  vary  in  height,  length,  and  breadth. 
Within  the  atom  area  they  are  known  as  wind  wares  and  often  appear  as  irregular 
mounds  (figure  1).  After  learing  the  storm  area  they  decrease  in  height  wd  inorease 
in  length  and  breadth  due  to  dispersion,  air  resistanoe,  risoosity,  turbulence,  eto. 
They  transferal  into  swell,  and  eventually  they  may  enter  shoaling  water  along  ooasts, 
finally  breaking  over  reefs,  against  cliffs,  and  along  beaohes. 

If,  during  an  amphibious  landing,  the  variation  in  breaker  heights  could  be  pre- 
dicted, a mechanism  would  be  available  by  whioh  landing  oraft  could  be  directed 
through  the  surf  during  intervals  of  relatively  low  breakers.  It  was  thought  that 
a possible  method  for  the  prediotion  of  relatively  calm  surf  might  be  formulated 
from  a study  of  the  relation  between  off-ehore  swells  and  the  breakers  on  a boaoh. 

If  suoh  a relationship  oould  be  eitabliehed,  readings  from  wave  reoord ere  operating 
aboard  off-shore  ships  could  be  used  during  days  ef  high  surf  to  direot  landing 
oraft  so  that  their  arrival  time  would  ooinoide  with  the  intermittent  intervals  of 
relatively  low  surf* 

Preliminary  studies  indicated  that  the  general  problem  was  extremely  complex  due 
to  the  dispersive,  short- ores ted,  and  refractive  properties  of  waves.  These  prop- 
erties are  discus aed  briefly  below i 

Diapers  Ion » Wave  velocity  in  deep  water  depends  primarily  upon  wave  length,  and 
to  a small  extent  upen  wave  height;  in  shoaling  water  wave  velocity  depends  upon 
wave  length  and  water  depth,  as  well  as  to  a small  extent  upen  wave  height.  The 
more  shallow  the  water  the  greater  is  the  effeot  of  depth  and  the  less  is  the  effeot 
of  wave  length  upon  wave  velocity.  As  long  as  the  velocity  is  dependent  upen  wave 
length  it  is  aaid  to  be  diaperaive. 


Consider  a typioal  reoord  suoh  as  the  one  shown  in  figure  2.  It  oan  be  seen  that 
both  the  heights  and  the  intervale  between  successive  orests  (hence,  lengths)  vary. 
Beoauee  ef  this  variability  it  oan  be  shown  in  a rather  simple  manner  that  the 
*vwavss"  oannot  be  of  permanent  form,  fbr  if  they  were,  then  the  orests  would  be 
moving  with  velocities  associated  with  the  lengths  and  ae  the  suooeeeive  leogtns 
vary,  the  velocities  would  vary;  henoe,  the  spaoes  between  oreste  would  ohange* 

It  follow*,  then,  that  they  oannot  be  of  permanent  form*  la  faot,  observations  ef 
waves  lb  the  generating  area  have  indicated  that  the  orests  oannot  be  followed  for 
more  than  a few  wave  lengths  i old  orests  gradually  disappear  cad-new  /ereetepv 
appear  (See  figure  7).  It  would  seem  reasonable  that  the  piateaoe . through' -whi«»h 
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individual  orosts  could  be  followed  would  dopoad  upon  the  width  of  tho  frequenoy 
spectrum.  Hence,  the  flrat  long  ewulls  arriving  from  a distant  storm,  consisting 
of  waves  of  nsarly  the  bias  frequecoy,  probably  oould  ba  traoked  for  a oonsidarably 
groator  distance  thin  the  wares  id  the  generating  area  (an  observed  faot).  The 
average  oondition  would,  of  course,  lie  between  these  two  limits. 

In  the  general  oaso,  then,  instead  of  following  individual  orests,  or  groups,  from 
one  looatlcn  to  another  and  merely  determining  the  travel  time,  one  must  treat  the 
transformation  in  both  spaoe  e-.i  time  of  a so-called  "oolored  noise*  random  time 
seriee.  In  order  to  deve1  ,p  the  relationship  between  the  disturbanoe  of  the  ooean 
eurfaoc  as  l;t  exists  offshore  and  the  breakers  along  a beaoh  it  is  neoeeaary  to 
utilise  Ibwier  Series,  Fourier  Integrals  or  random  funotion  analyses  of  the  surface 
fluctuation,  combined  with  the  linear  wave  theory. 

Fbrtunately,  the  dependenoy  of  wave  velooity  on  wave  length  decreases  as  the  water 
depth  deoreMes.  In  faot,  on  a beaoh  the  velooity  of  a wave  as  it  breaks  is  nearly 
independent  of  lengthi  it  depends  almost  entirely  upen  the  water  depth.  Thus, 
an  observer  on  a pier  watehing  relatively  long  swells  travel  from  the  end  ef  the 
pier  to  tho  breaker  sons,  can  easily  follow  the  individual  oreets.  Because  of  this 
relative  insensitivity  near  shore  it  is  felt  that  a relatively  simple  solution  may 
be  obtainable  for  many  oases. 

Short-oreetedneee  (See  Figure  8)1  Consider  the  oase  of  a wave  reoorder  being  located 
as  is  shown  in  Figure  4.  Because  the  waves  are  short-oreeted,  the  height  of  the 
breakers  oan  be  dlrsotly  predioted  only  along  a relatively  short  seotion  of  the 
beaoh*  It  ether  leoertions  the  breaker  heights  may  be  higher  or  lower.  Fbr  beaohes 
of  snail  length  this  will  not  be  ae  important  as  for  long  beaohss. 

Observations  indioate  that  the  oreet-length  (i.e.,  breadth)  depends  largely  upon 
the  dietanoe  the  waves  have  travelled  from  the  storm  area;  the  greater  the  distance 
the  longer  the  oreet-length,  Eenoe,  swell  from  distant  storms  should  present  less 
a problem  than  looal  waves. 

Refraction;  Ihen  waves  enter  water  in  wkioh  the  depth  is  less  than  half  their 
length  't'&eir  velooity  depends  upon  the  depth  as  well  as  their  length,  so  waves 
travelling  at  an  angle  to  the  underwater  contours  will  bend  (Figure  4).  Beoause 
of  this  it  is  not  possible  to  prediot  dlrsotly  the  offset  of  wave  transformation 
along  a lino  parallel  to  their  direotion  ef  travel  in  deep  water;  refraotioa 
must  be  considered. 

Actually,  the  problem  ie  even  more  dlffleult  beoause  of  the  dispersive  property 
of  waves,  oembined  with  their  non-uniform  oharaoterlstios.  The  various  components 
refract  in  sua  appropriate  manner;  henee  the  wave  motion  on  the  water  surf&oe  le 
transforming  in  a direotion  perpsndioular  to  the  direotion  of  propagation  as  well 
as  in  the  direotion  of  propagation.  Fortunately,  in  the  more  shallow  water  near 
beaohes  this  dispersive  problem  beoomes  less  aoute. 

Considerable  theoretical  studies  have  been  made  in  regard  to  the  preblems  due  to 
the  dispersive  properties  of  ooean  waves.  These  studies  have  been  compared  with 
results  obtained  in  the  laboratory.  In  addition,  arrangements  wsrt  made  with 
the  Santa  Crus  Portland  Cmeat  Company  for  ths  use  of  their  pier  (looated  at 
Davenpert,  California)  - Figure  12,  and  four  step-res  is tant  gages  wsre  ordered 
from  the  Beach  Erosion  Board  for  installation  along  ths  pier.  These  have  Just 
beta  delivsrtd,  and  a detailed  field  program  will  oommenoe  shortly. 
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Figure  2 


COPY  OF  WAVE  RECORD  FROM  MARK  2 PRESSURE  HEAD 
IN  20  FEET  OF  WATER,  SOLDIERS  CLUB  BEACH,  FORT  ORD,  CALIFORNIA 


Be  LABORATORY  STUDIES 


gave  Channel  > The  laboratory  experiments  war®  performed  la  a waT®  channel,  1 foot 
wide,  60  f«st  long  and  3 faot  deep  looated  la  the  Fluid  Meohanloa  Laboratory  of 
the  University  of  California,  Berkeley*  The  front  side  of  the  ohaanei  consisted 
of  plate glass,  framed  in  3 ft*  x 3 ft.  steel  frames,  through  which  motion  pictures 
of  the  wave  transformation  were  takon.  The  wave  generator,  of  the  flapper  type, 
was  looated  at  oae  end  ef  tb»  onannel*  Both  the  amplitude  and  period  of  the  flapper 
movement  were  adjustable.  The  period  of  the  flapper  movement  oould  be  varied  be- 
tween approximately  0.4  seconds  and  2 seooads.  At  the  oppesite  end  of  the  ohannel 
from  the  wave  generator  an  aluminum  beaoh  was  installed  for  the  purpose  of  elimin- 
ating refleoted  waves* 

Transformation  of  Finite  Have  Groups  in  a Uniform  Dspth  of  later t lave  groups 
were  generated  by  operating  the  flapper  manually  for  the  shorter  wave  groups  and 
meehanioally  for  the  longer  groups*  Manual  operation  gave  both  better  oontrol  of 
the  number  of  waves,  and  a more  uniform  wave  period.  It  was  found  when  operating 
the  generator  meehanioally  that  it  took  approximately  B waves  before  the  generator 
was  aooeleratod  to  a constant  period,  and  that  at  the  end  of  the  run  the  flapper 
did  not  stop  ae  soon  aa  the  motor  was  turned  off.  Hence,  meohanloal  operation  of 
the  flapper  resulted  in  wave  groups  with  long  waves  at  the  start,  then  & gradual 
deorease  in  wave  length  (beooming  oonatant  in  lungth),  and  then  after  the  motor 
wae  stopped,  a gradual  inorease  in  wave  length  again,  until  the  flapper  completely 
stopped* 

The  initial  groups  studied  contained  from  1 to  25  waves.  The  surfaoe-time  histories 
were  recorded  at  4 locations  along  the  channel  with  the  first  element  1 foot  from 
the  generator,  the  seoond  11  feet,  the  third  21  feet  and  the  fourth  31  feet.  The 
water  depth  was  uniform  and  equal  to  2.0  feet  (at  elements  1,  2 and  3)  from  the 
flapper  to  the  toe  of  the  beaoh.  It  was  1.6  feet  at  element  4,  which  was  looated 
landward  of  the  toe  of  the  beaoh.  The  beaoh  slope  was  1 s 15. 

The  aurfaoe-time  histories  wers  recorded  by  use  of  double-wire  resistance  elements* 
connected  to  Brush  reoorders.  The  data  on  these  reoorders  were  correlated  by 
connecting  them  to  & switch  operated  by  the  wave  generator,  at  station  0+00.  The 
flapper  passing  the  switch  the  first  time  closed  it,  causing  a murk  to  be  placed 
on  the  reoorda,  thus  affording  a "aero  time".  Figurs  5 illustrates  one  of  the 
runs  (Ho*  94)  with  the  flapper  of  the  wave  generator  moved  twice  back  and  forth, 
with  a result  of  two  waves  is  the  initial  group  one  foot  from  the  generator.  In 
this  figure  the  transformation  of  the  wave  group  as  it  advanoes  down  the  ohannel 
can  readily  be  seen* 

Knowing  the  surface  elevation-time  history  at  one  location  it  is  possible  to  pre- 
dict the  water  surfaoe  at  the  other  locations.  This  has  been  done  and  the  results 
oompared  with  the  theory,  with  good  agreement.  The  results  are  presented  in  sub- 
sequent sections  of  this  report. 

Wave  Transformation  of  Son-uniform  Wavee  in  Uniformly  Shoaling  Water;  The  experi- 
mentis  just  described  were  performed  to  check  the  theory  for  prediction  in  watey 
of  uniform  depth.  But  the  important  problem  in  most  practical  canes  is  to  prediot 
wave  trana format ion  in  shoaling  watsr.  Hence,  a beach  with  a uniform  slope  of  lg40 

*Ho?ieon,  J«R.  "Measurements  of  Heights  by  Resistance  Elements”.  The  bulletin, 
Beaoh  Erosion  Board,  Corps  of  Engineers,  vol*  3,  no.  3,  pp.  16-22 ? i94S* 


B. 


m • ina tailed.  Beoauee  of  the  limited  leagth  of  channel  a watar  dapth  of  only  1,17 
faat  was  uaad.  This  la  ft  a seotion  of  uniform  dapth  10.42  feat  long  between  the 
toe  of  the  beaoh  and  the  ware  generator,  which  allowed  the  wares  to  stabilise  be- 
fore reaching  the  toe  of  the  boaoh.  The  water  line  on  the  beach  was  57.80  feet 
from  the  ware  generator, 

WaTa  groups  were  generated  by  moving  the  flapper  manually  in  order  to  mi*  periods 
and  amplitudes.  In  addition  to  the  finite  ware  groups,  measurements  were  made  for 
steady-state  uniform  wares  by  ^aerating  the  wares  meohanioally  and  taking  the 
measurements  only  after  a - ieady  state  was  reached,  The  water  surface-time  history 
was  reoord^d  at  two  lections.  The  first  element  was  looated  just  landward  from 
the  toe  of  the  beaoh  (10,42  faat  from  the  flapper)  with  the  water  depth  d^  - 1.17 
feet.  The  seoond  element  was  looated  41,62  feet  from  the  flapper,  with  d - 1/3  - 

0,39  foot. 

As  far  ae  the  uniform  trains  of  wares  were  aonaarned,  it  was  easy  to  identify  the 
same  ware  on  both  reoords  and  to  determine  the  number  of  wares  between  the  elements. 
However,  diffioultles  arose  in  the  ease  of  traino  of  non-uniform  waves.  Some  of 
the  waves  would  disappear  and  new  waves  appear;  so,  it  was  not  always  possible  to 
identify  the  same  wave  on  different  reoords.  Unfortunately,  reoords  of  this  sort 
do  not  show  the  continuous  ohange  of  a train  of  irregular  waves  even  when  the  points 
of  measurement  are  numerous  and  oloee  together.  As  this  was  very  important  to  the 
understanding  of  the  problem,  it  was  dsoided  to  obtain  a oonMuuous  reoord  of 
irregular  waves  with  movies.  Pour  35  mm.  Sell  4 Howell  "Bymo*  cameras  (25  mm. 
fooal  length  f.  2.S  lenr)  were  used,  eaoh  one  having  a coverage  of  9 feet.  The 
cameras  were  operated  simultaneously,  overlapping  eaoh  other's  field  of  view. 

In  order  to  obtain  a time-aoale  for  the  measurements t and  to  oorr elate  the  data 
of  the  different  oameras,  three  olooks  were  used.  Two  of  the  olooks  were  gradu- 
ated in  0.01  seoend  increments  (one  sweep  of  the  arm  equal  to  1 seoond)  and  one  in 
0,01  minute  (l  sweep  equal  to  1 minute).  The  arrangement  of  the  set-up  is  shown 
in  Plgure  6.  The  olook  graduated  in  0,01  minute  increments  (clock  2 in  the  sketch) 
was  used  only  to  oorrelate  t.ae  readings  of  olooks  1 and  3. 

The  measurements  obtained  fVx&  the  motion  pictures  were  oheoked  by  placing  a 
double  wire  resistanoe  el  am  out  in  the  view  cf  Camera  III  at  a known  location,,  and 
measuring  the  surfaoe-time  h:. story  by  means  of  a Brush  recorder  at  the  same  time 
that  the  movies  were  being  Is. ken.  To  oorrelate  the  data  of  the  Brush  recorder 
and  the  movies,  the  oonnaotiens  to  the  olooks  were  made  through  the  Brush  recorder 
so  that  the  moment  the  olooks  were  started  a mark  was  plao&d  on  the  record.  Agree- 
ment between  the  data  obtained  from  the  movies  and  that  from  the  Brush  reoord  waa 
found  to  be  very  good. 

Experimentation  with  the  photography  showed  that  clear  water  gave  a surface  line 
in  the  photographs  whloh  was  difficult  to  read.  There  appeared  to  be  several 
methods  of  overoeaniag  this  and  thereby  obtain  a clear  image  of  the  surface  profile 
in  the  piotures.  It  was  decided  to  oover  the  glass  windows  in  the  channel  with 
tracing  paper,  stretched  tightly  against  the  surfaoe  of  the  glass.  A grid,  with 
OoCl  foot  divisions,  was  drawn  on  the  paper  to  obtain  a scale  for  the  evaluation 
cf  the  data,  and  the  distunoee  from  the  wave- generator  indicated.  When  a strong 
light  was  directed  to  the  water  surfaoe  a vory  clearly  distinguished  shadow  line 
was  obtained  on  the  tracing  paper;  henoe,  the  shadow  of  the  water-surfaee  profile 
was  aotually  photographed.  The  results  w era  satisfactory.  Special  care  was  taken 
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to  keep  ths  traoing  paper  aa  tight  as  pvuaible  to  the  glass;  otherwise  the  shadow 
image  would  not  be  dear  and  would  result  in  erroneous  readings.* 

The  duration  of  th*  run  was  approximately  90  seconds , as  this  was  oonsidered  to 
b«  adequate  from  a statistical  standpoint. 

The  movies  wore  analysed  frame  by  frame,  eaoh  war®  in  the  photographs  being  assigned 
a number,  and  the  time  and  location  of  these  waves  being  measured  as  the  waves  ad- 
vanced down  the  ohannel.  The  prooedure  was  very  tedious  and  the  evaluation  of  the 
data  has  not  yat  been  conflated.  The  data  obtained  were  plotted  as  time  versus 
distance  traveled.  H^oe,  the  diatanoe  between  waves  in  the  vertioal  direotion 
gives  'the  wave  period*  at  any  desired  time  and  looation,  the  distanoe  between  waves 
in  the  horisontal  direotion  gives  the  wave  lengths,  and  the  slopes  of  the  ourves 
give  the  wave  velocities,  A sample  of  the  data,  from  Camera  I,  is  shown  in  Figure 
7.  These  are  the  waves  in  the  deepest  water;  henoe  the  tlme-distanoe  ourves  are 
almost  straight  lines.  It  is  expeotsd  that  the  curvatures  would  increase  as  they 
oome  into  shallower  water  and  the  velooities  of  the  wave  travel  decrease. 

The  evaluation  of  the  data  is  not  yet  advanoed  enough  to  disousn  fully  the  resultB, 
but  looking  at  FLgure  7 it  is  interesting  to  note  that  when  a wave  disappears,  there 
appears  a new  wave  a short  diatanoe  to  the  rear  of  the  previous  wavs.  This  is 
demonstrated  by  Wave  88,  whioh  disappear*  at  Station  33  and  is  replaced  at  8tation 
34  by  Wave  87-A,  The  same  oharaoteristios  oan  be  found  for  Waves  87  and  83.  Wave 
87  breaks  and  disappears  and  ia  replaoed  a short  distanoe  later  by  another  wave. 

Waves  88  and  87-A,  and  87  and  86  demonstrate  the  mixing  of  two  waves.  In  the  aaeo 
of  Wares  88  and  87-A,  the  two  oresta  merge,  beooming  a long,  flat-orested  wave, 
with  no  distinguishing  maximum  (see  Figure  7,  Stations  33-34).  The  single  crest 
then  separates  into  two  oresta  again, one  of  whioh  soon  disappears.  In  the  oaaa 
of  Waves  87  and  86  a different  phenomenon  oooura  i Wave  87  breaks,  aid  water  rolls 
forward  merging  with  Wave  86  and  forming  a long,  flat-oreBted  wave.  As  in  the 
oaoe  of  Waves  88  and  87-A  this  separates  into  two  waves,  one  of  which  soon  disappears. 

In  Figures  8-10  are  shown  a series  of  skstohes  of  the  wave  profiles  as  they  moved 
dawn  the  ohannel.  Eaoh  aketoh  shows  nearly  the  entire  shoaling  seotion.  The  time 
interval  between  eaoh  sketoh  varies. 

Wind  Waves  in  Uniformly  Shoaling  Water i Several  laboratory  experiments  were  per= 
formed  in  order  to  determine  the  oharaoteristios  of  wind  generated  waves.  The 
ohannel  was  covered  to  fora  a wind  tunnel  with  the  water  as  the  bottom  boundary 
of  the  wind.  A sloping  beaoh  was  plaoed  at  one  end  of  the  channel.  The  wind  was 
foroed  through  the  tunnel  over  the  water  and  toward  the  beaoh.  The  resulting  wares 
were  measured  simultaneously  at  three  stations  along  the  sloping  beaoh.  Figure  11 
shows  a sketch  of  the  experimental  set-up. 

The  analyaiB  of  these  data  will  show  the  effect  of  shoaling  on  wind  waves  together 
with  their  transformation  in  the  generating  area  due  to  their  dispersive  properties. 
The  results  will  he  presented  in  a subsequent  report. 
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WAVE  TRANSFORMATION  - IRREGULAR  WAVES 
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SCHEMATIC  DIAGRAM  OF  EXPERIMENT  ON  WIND  WAVES  IN 
UNIFORMLY  SHOALING  WATER 


C.  FIELD  STUDIES 
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Fit  Id  Work  at  Davenport,  California!  Through  the  oooperation  of  the  Santa  Crus 
PortlacZ  Semen?'  Company,’  the  University  of  California  has  been  given  permission 
to  use  an  abandoned  oement  loading  pier  at  Davenport,  California  for  studying 
wave  transformation  (Figure  12).  The  pier  is  unusually  long  (2200  feet)  and  ie 
oonstruatad  in  an  exposed  looatlon  on  the  ooaat,  giving  an  excellent  ohanoe  for 
the  study  of  long  period  large  vvaves  from  the  North  Paoifio  area. 

Moving  ploturea  were  ta*.en  of  waves  moving  along  the  pier  in  an  initial  effort  to 
define  the  oomplexity  of  the  transformation  problem.  It  was  found  that  waves 
oould  be  traoked  visually  from  the  end  of  the  pier  to  shore  and  that  the  travel 
times  of  these  wuves  oould  be  predioted  satisfactorily  from  theoretlo&l  consid- 
erations. Work  is  now  going  on  in  preparation  for  an  extensive  wave  measurement 
program  correlating  waves  along  the  pier  and  waves  a mile  westward  of  the  pier  end. 
It  is  planned  to  install  four  Beaoh  Erosion  Board  step-re sistanoe  gages  along  the 
pier  to  reocvd  the  transformation  of  waves  in  shallow  water  and  to  install  a net- 
work of  wave  pressure  reoorders  one  mile  off  shore  is  a depth  of  80  feet  to  measure 
the  deep-water  waves  ooming  inco  shore. 

It  has  been  uooessary  to  install  naw  walkways  along  the  pier  and  also  to  put  in 
eleotrio  service  and  oonmunioatior  lines  sinoe  some  of  the  pier  supers truoture 
has  rotted  or  washed  away  and  muoh  of  the  original  equipment  has  been  removed. 

A wave  pressure  reoordar  was  put  into  operation  at  the  pier  end  to  give  some 
idea  of  the  wave*  site  and  period  to  be  ezpeoted  during  the  proposed  tests. 

Field  work  at  Bllwoodt  In  oonneotion  with  the  problem  of  wave  transformation, 
two  preliminary  tests  were  made  at  Ellwood,  California  on  & pier  owned  by  Si  goal 
Oil  and  Oas  Company.  The  tests  oonsieted  of  spacing  two  wave  pressure  reoorders 
a distance  of  28. b and  205  feet  apart,  near  the  surface  and  on  & line  with  the 
direotion  of  approaoh  of  the  waves.  The  water  depth  at  the  ueaward  reoordar  2200 
feet  from  shore  was  36  feet.  Analyses  of  the  data  obtained  from  these  tests  are 
deaoribed  in  Sections  F and  0, 
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D.  Lima  PREDICTION  IN  WATER  OP  CONSTANT  DEPTH 


In  a previous  report*  the  predioted  surface  elevation  at  a diatanoe  x in  the 
direction  of  wave  propagation  was  shown  to  be  of  the  form 
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where  K 
deep  water  we  find  explicitly 

X'.t)  = - 


which  is  plotted  vs.  t for  x - 10  In  Figure  13*  In  general  this  surfaoe  kernel 
o&n  be  written  in  the  fora 
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Calculations  have  been  performed  for  4 : 6 and  y®  : 2,w©  . This  last  set  of 

2ir  rr 

values  implies  the  neglect  of  component  periods  of  less  than  yj  — or  the 

consideration  of  d/L  values  less  than  0,64, oo  respectively.  The  integral  has 
been  calculated  using  Filoa's  formula.  The  results  are  plotted  in  Figure  13. 


Suooessive  predictions  were  made  for  d a 2 feat  and  stations  10  feet  apart  the 
first  being  one  foot  from  the  wave  paddle.  The  results  are  indioated  in  Figure 
14  for  an  Initial  finite  sine  wave  the  first  record  being  taken  one  foot 
from  the  wave  flapper.  The  results  at  Station  u were  predicted  from  the  measured 
surfaoe-time  history  at  Station  2 and  the  prediction  for  Station  4 was  made  from 
the  measured  surfaoe-time  history  at  Station  3,  In  Figure  15  a single  prediction 
is  made  for  a larger  number  of  initial  oscillations.  Trial  predictions  were  made 
for  data  secured  at  Ellwood  but  since  the  experimental  conditions  did  not  con- 
form to  those  used  in  the  calculations  no  comparison  will  be  offered. 
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S.  LXJJlkAR  FBBDICTXON  OP  WAVS  IRA  N 3 FORMAT  I QH  OK  BSACHES  OF  COKSXANT  SLCjPg 


In  order  to  generalise  the  prediotion  of  ware  transformation  from  water  of  oon- 
stant  depth  to  ahoaling  water  one  must  in  general  take  aooount  of  ohangea  in 
phase,  amplitude  and  refraotion  angle.  Fbr  simplicity  we  shall  oonaider  beaohes 
of  uniformly  oonetant  olope  with  wawes  mowing  direotly  onahore  and  hawing  areata 
parallel  to  the  bottom  oontoura.  We  aaaume  that  a periodio  ware  has  an  ampli- 
tude whioh  ia  determined  by  Rayleigh's  hypotheaia  of  the  oonaorwation  of  trans- 
mitted power  and  the  phaa*  ia  determined  by  suitably  Integrating  the  instantaneoua 
wawe  number  appropriate  to  the  depth.* 


we  find 
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Pbr  a periodio  wawe  of  frequeaoy  the  surfaoe  profile  is  of  the  form 
*f\  (t)  = a (t)  ooi  ( \ (d)  - O’  t) 

where 
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x being  measured  plus  toward  shore  with  x a o at  depth  dg. 


24. 


For  a continuous  sequenoe  of  underlying  frequencies  we  generalise  thia  simple 
formula  by  integration  over  C in  order  to  gain  Fourier  integral  representations, 
Thus 


^2  (t)  = j a2  (S')  d«- 
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where  IT  " - g k tanh  V u,  Inverting  (1)  by  Fourier's  integral  theorem  the 
spectral  amplitude  is 
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Applying  the  method  of  stationary  phase  to  this  integral  results  in  tho 
approximate  formula 
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Calculations  using  this  approximation  were  performed  for  - l/40,  dx  - 0,39  ft. 
and  do  - 1.17  ft.  corresponding  to  experimental  conditions  described  elsewhere 
in  this  report.  Numerical  oaloulatioas  were  oarried  out  for  the  range  of 
periods  between  0.36  and  1.26  seo.  which  appeared  to  oover  the  range  of  experi- 
mental  conditions.  The  results  are  given  numerically  in  Table  1,  The  amplitude 
f (r)  and  the  phase  g(vj  are  plotted  separately  in  Figure  16.  Predictions  employ- 
ing this  kernel  are  compared  with  experiment  in  Figures  17a  and  b.  Discrepancies 
between  theory  and  experiment  ar6  in  part  attributable  to  the  approximations 
involved  in  employing  the  method  of  stationary  phase. 
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If  the  linear  theory  for  wave  propagation  ie  adequate  it  should  be  possible  to 
start  with  a record  of  the  surface  time  history  of  wave  motion  at  some  point, 
and  then  after  harmonically  analyzing  a section  of  this  recor'  into  its  various 
components,  propagating  the  components  separately  in  the  direction  of  travel 
to  some  desired  point  and  then  recombine  them.  The  point  of  reference  on  the 
longest  component  would  arrive  first  at  the  point  of  interest  while  the  ref- 
erence point  of  the  shortest  component  would  arrive  last.  The  recombined  surf- 
ace time  history  would  begin  -o  be  valid  only  when  the  shortest  components  had 
started  to  arrive  at  the  point  of  interest?  hence  the  interval  for  which  the 
surface  time  history  can  be  predicted  at  this  second  point  is  always  less  than 
the  length  of  record  at  the  first  point.  If  the  section  of  record  at  the  first 
point  is  relatively  short  and  the  distance  between  the  two  pointe  is  relatively 
long  then  no  prediction  from  the  first  to  the  second  is  possible.  This  indicates 
that  "■hs  greater  the  distance  to  the  point  of  prediction  the  greater  would  be 
the  section  of  measured  record  that  must  be  harmonically  analysed.  It  follows 
that  the  longer  the  section  of  record  harmonically  analyzed  the  greater  the 
number  of  harmonics  necessary. 

The  basic  limitations  for  the  harmonic  analysis  method,  as  presented  herein, 
are  that  the  depth  of  water  stays  constant  (the  effect  of  shoaling  is  not 
included)  and  that  the  waves  behave  according  to  the  linear  wave  theory.  The 
mathematics  involved  in  this  method  start  with  the  representation  of  the 
measured  surface  by  a Fourier  Series  which  is  given  by  the  expressions 
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where 


??]_(t)  = measured  or  known  surface  time  history  at  point  1 

”7 2 (t)  * computed  surface  time  history  at  point  2 or  any  point  other 
than  point  1 

aQ,  a^,  ...  an,  bp  ....  b - Fourier  coefficient 
t - time  - seconds 

n - number  of  sine  or  cosine  components  which  is  one-half  the  number  of 
ordinates  measured  in  the  type  of  harmoii-LC  anaysis  used* 

T s length  of  record  harmonically  analyzed  in  seconds 

IS  ■ 'T ; T«  « -jr  .....  T = -l—  = w period  of  the  component  in 

22  n n 8300nd8 

A = distance  between  point  1 and  point  2 or  any  other  point  measured  in 
direction  of  wave  travel.  Points  must  be  in  line  in  this  direction. 


*Den  Sartofy-J.P.  "Mechanical  Vibrations",  McOraw-Hill  Book  Co.,  Inc. 

New  York  and  London,  191)0 
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where 


= t - t]_  time  scale  on  record  2 in  seconds 
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"jT*  5 t2  = “q~  J ...  tn  = -4j-  time  for  each  component  to 

*2  ^n  propagate  the  distance  ^ 


°1  = Tx  ; 


- wave  velocity  of  each  component 


Lx  = wave  length  of  each  component  obtained  from  wave  theory  for  known 
d and  Tn* 

d = still  water  depth  ft. 

Figure  18  shows  a sketch  of  the  time  scale  used  with  regard  to  (t'+tx-tn). 

At  position  2 the  time  of  each  component  has  a definite  angle  or  argument 
g'fft 

(“t — ) which  for  evaluation  purpose  must  be  separated  into  the  time  argument 
xn 

(—m—  ) plus  a phase  angle  This  later  combination  was  used  as 

xn  rn 

^oriofigi|^U^>^^Mlysis  of  Subsurface  Pressure  Records  in  Constant  Depths  and  on 
Slopjjjqfrjtl—w — 1 *y»8eries  3,  Issue  336,  Institute  of  Engineering  Research,  Univ- 
ersity of  California,  Berkeley.  California,  May  1952.  and  Wiegel,  R.L.,  "Oscillatory 
Waves'1,  Bulletin  of  Beach  Erosion  Board,  Spec.  Issue  1,  19uo,  Corps  of  Engineers 
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describe  1 in  Equation  2.  The  resulting  Equation  3 may  be  simplified  mathematically, 
but  such  action  may  tend  to  complicate  the  final  application  which  is  relatively 
simple. 

Several  experiments  were  conducted  in  both  the  laboratory  and  the  field.  A few 
of  these  were  analysed  with  various  degrees  of  success.  In  order  to  obtain  good 
results  the  experiments  and  application  of  the  method  must  follow  or  adhere  to 
several  rules  as  follows i 

(i)  The  points  zi  measurements  and  prediction  must  be  in  line  and  in  the 
direction  of  wave  propagation. 

(ii)  The  wave  system  must  be  regular  enough  to  have  only  one  direction  of 
propagation.  The  success  of  the  method  depends  upon  the  waves  behaving  according 
to  the  linear  wave  theory. 

(iil)  At  least  three  to  six  ordinates  should  be  measured  on  the  longest  wave 
in  order  to  obtain  a satisfactory  harmonic  analysis.  Since  only  the  1*8  point 
harmonic  analysis  aoheme*  was  uaad  this  would  limit  the  length  of  record  analysed, 
the  greatest  distance  to  the  predicted  point  and  the  length  of  the  predicted 
record. 

There  are  of  course  other  ways  of  doing  this  type  of  transformation  but  inherent 
to  all  is  the  limitation  that  the  greater  the  distance  propagated  the  greater  is 
the  length  of  measured  record  that  must  be  analysed.  It  might  be  better  to  say 
that  this  limitation  is  at  least  inherent  to  the  linear  wave  theory. 

Figure  19  shows  the  results  of  the  laboratory  study  of  wave  transformation. 

Figure  19b  shows  the  original  measured  surface  wave  system  and  its  reconstruction 
by  harmonic  analysis  which  is  a check  on  the  accuracy  of  the  wave  components. 

Figures  19c,  19d,  and  19e  show  the  transformation  of  the  wave  after  a travel 
distance  of  10,  20  and  30  feet  respectively  together  with  the  predicted  or  com- 
puted wave  time  history  by  the  harmonic  analysis  method  or  Fourier  Series  analyBis. 

Figure  20  shows  a field  study  where  Figure  20b  is  the  measured  wave  pressure 
and  its  reconstruction  by  harmonic  analysis  as  a check  on  the  wave  components. 

Figure  20c  shows  the  wave  transformation  after  201*  feet  of  travel  together  with 
the  predicted  wave  transfo nation 

The  results  shown  here  are  fairly  good,,  but  not  all  attempts  were  this  accurate. 
Additional  analyses  of  the  accurately  controlled  field  tests  in  progress  will 
help  to  determine  the  validity  of  this  method.  One  indication  at  this  time  is 
that  the  greater  the  travel  distance  the  greater  will  be  the  error.  This  error 
appears  in  the  phase  relationship  of  the  various  components.  Of  course  the 
effect  of  shoaling  has  not  been  taken  into  account.  This  problem  is  mainly 
whether  or  not  all  components  can  be  separately  transformed  by  shoaling  and 
then  recombined  to  give  the  actual  transformed  wave*  This  of  course  implies 
that  the  linear  wave  theory  holds  for  shoaling  waves,  which  is  not  true  in 
general  in  relatively  shallow  water. 


*Den  Hartog,  J.P.  "Mechanical  Vibrations",  McGraw-Hill  Book  Co,,  Inc., 
New  York  and  London,  19 1*0 
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Introduction 


Within  the  problem  of  wave  transformation  there  may  arise  the  question  of  the 
degree  of  predictability  inherent  in  the  phenomenon  apart  from  the  theory 
employed  for  making  the  prediction  or  the  means  of  carrying  it  out.  If  in 
particular  a finite  linear  '..-mbination  of  discrete  equally-spaced  values  of 
one  observed  variable  ir  to  be  used  to  predict  the  other  variable,  the  method 
of  least  squares  may  bo  employed  to  yield  not  only  an  optimum  prediction  kernel 
of  this  type  but  also  a measure  of  the  associated  error  of  prediction.  The 
preliminary  studies  described  below  were  undertaken  in  order  to  show  the  de- 
gree of  accuracy  possible  in  predicting  the  transformation  of  waves. 


The  least  squares  method  will  rest  on  the  interpretation  of  the  two  time- 
hi8toriee  f^*)  as  the  components  of  a two-dimensional  discrete  station- 

ary random  process.  That  is,  they  constitute  a realization  of  a probabalistic 
system  characterized  (for  present  purposes)  by  three  functions  yC-*-*), 

V/(22)j  aach 

being  the  correlation  function  associated  with  fv^)  and 
The  direct  computation  of  the  least-squares  prediction  kernel  described 
below  depends  on  the  previous  knowledge  of  certain  ordinates  of  the  correlation 
functions  and,  while  conceptually  simple,  its  exact  computation  becomes  tedious 
when  the  length  of  the  kernel  exceeds  more  than  a few  terms.  The  method,  how- 
ever, is  applicable  to  any  two  time-history  functions  provided  only  that 
together  they  may  be  regarded  as  a random  process  of  the  kind  described. 

The  required  values  of  the  correlation  functions  must  be  estimated  from  the  ob- 
served sample  time  history,  it  being  assumed  that  this  sample  represents  a 
certain  infinite  (unobservable)  time  history. 


The  length,  n,  of  the  prediction  kernel  employed  in  the  preeent  studies  was 
dictated  by  computational  convenience.  The  time  lag,  l , between  the  time 
of  the  predicted  ordinate  f(*)  and  that  of  the  latest  ordinate  f^1'  used  from 

the  observed  time  hiatory  was  determined  so  as  to  yield  the  best  prediction. 
The  prediction  kernel  then  consists  of  a sequence  of  numbers  dQ,  do.  ...,  dn_]_ 
which,  when  applied  to  the  observed  time  history,  produce  dk  fW  as  the 

value  predicted  for  f(2)  . mo 

tti 

Linear  Least-Squares  Theory 


It  is  assumed  that  suitable  constants  have  been  subtracted  from  the  values  of 
each  of  and  f(2)  so  as  to  reduce  the  mean  value  of  each  to  zero.  Being 
given  the  prediction  lag  & and  the  kernel  length  n,  the  problem  is  to  determine 
the  sequence  of  real  number’s  d0,  d]_,  ...,  dn_]_,  such  that  when  the  linear  combin- 


ation 


is  formed  for  each  given  value  of  t,  the  mean,  or  expected, 
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value  of  the  squared  error  $2  = _ jijmJ  f^^^J 
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for  j * o,  1,  . ..,  n-1,  the  bar 


denoting  the  (linear)  mean-value  operator.  Introducing  the  correlation  functions 
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this  has  the  form 
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n-1. 


The  conditions  on  the  d^  are  thus  those  imposed  by  q simultaneous  linear 
equations  in  n unknowns.  Assuming  these  have  a unique  solution,  the  corres- 
ponding prediction  kerne"1  *111  yield  the  minimum  least-squares  error  given  by 


min 
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ueJng  the  condition  on  the  'a.  A relative 


measure  of  the  error  of  prediction  suitable  for  many  purposes  1b  provided  by 
the  ratio  of  the  root-mean-aquare  error  to  the  root-mean-square  deviation  of 
the  observed  time  history  from  its  average  value.  This  ratio  may  be  written 
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and  may  be  computed  from  the  values  of  the  correlation  functions  alone,  It 
is  thus  possible  to  determine  what  prediction  lag  and  kernel  length  may  best 
be  used  before  any  actual  predictions  are  made. 


Applications  of  the  Theory 

The  finite  linear  least-squares  prediction  kernel  was  computed  for  three  ex- 
perimental seta  of  data,  each  corresponding ’to  the  transformation  of  wave 
systems  between  two  points  displaced  norizontally  in  the  general  direction  of 
wave  propagation. 


I.  laboratory  Wave  Channel  - Here  the  rave-generating  flapper  was  moved 
briefly  and  a short  wave  train  of  approximately  one-second  waves  was  started 
down  the  wave  channel.  At  four  stations,  each  one  approximately  ten  feet  from 
the  preceding  one,  time  histories  of  the  surface  elevation  of  the  water  were 
recorded  for  approximately  30  seconds.  Predictions  were  made  from  the  first 
station  to  each  of  the  other  three,  as  well  as  from  the  second  to  the  third. 

Two  representative  examples  are  described  below. 


A.  Station  Li  from  Station  1 - From  measurements  made  every  0,1 


second,  values  of 


(12) 

112 


y(lD 

0 o » 


yU!) 
a ]_  s 


yU2) 
« 106, 


were  estimated  and  used  to  compute  an  eight-term  prediction 


kernel  do,  d]_,  dj,  Here  X - 105,  corresponding  to  10,5  seconds, 

roughly  ten  times  the  wave  period.  The  resulting  prediction  may  be  com- 
pared wj.th  the  observed  values  at  the  second  point  in  Figure  21a. 


) 


*The  non-normalized  correlation  function  for  anv  two  time  history 

functions  fd)f(J)  j_8  defined  to  be  the  average  product  f^r^p.  It  Trill  be 

eeen  that,  for  i“J,  and  when  the  lag  difference  between  subscripts  is  zero,  the 
c^^yelation  function  is  merely  the  square  of  the  r.m.s.  value  of  the  function 
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Station  3 from  Station  2 - From  measurements  made  every  0.1  second, 

(11) 
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B. 


values  of 


y (ID 

0 o » 


y(ll) 

0 1 » 


and 


y/(12) 

0 21, 


y(12) 

o 22  > 


(12) 

28 


were  estimated  and  used  to  compute  an  eight- term  prediction  kernel  d0,  d^, 
d^.  Here  i * 21,  corresponding  to  2.1  seconds,  roughly  twice  the  wave 
period.  The  resulting  prediction  may  be  compared  with  the  observed  time 
history  at  the  second  point  in  Figure  21b. 


II.  Ocean  Waves,  28.5-foot  Separation  - Here  the  surface-elevation  time  history 
due  to  natural  ocean-wr-  d 'conditions  was  measured  for  a period  of  approximately 
five  minutes  at  two  stations  separated  by  28.5  feet.  Predictions  were  made  from 

measurements  made. every  0," 


the  seaward  station  to  the  shoreward  station.  From 
second,  values  of 

were  estimated  and  used  to  compute  a four-term  kernel  d0,  d^,  dg,  d-j 


and 


Here  & was  taken  as  2,  corresponding  to  1.0  second,  roughly  one-tenth  of  the 
mean  was  period.  The  resulting  prediction  may  be  compared  with  the  observed 
time  history  at  the  second  point  in  Figure  22. 


III.  Ocean  Waves.  205-foot  Separation  - Here  the  subsurface  fluctuation  time 
history  due  to  natural  wave  conditions  was  measured  for  a period  of  1000  seconds 
at  two  points  separated  by  205  feet  and  at  depths  of  approximately  35  feet.  Pre- 
diction was  made  from  the  seaward  station  to  ^xshorey^  stat^cj^  Froip^as- 


urements  made  every  0.25  second  values  of 


in- 


and  ^23^»  ^^35  were  e8tinated  ^ used  to  compute  the 

four-tern  prediction  kemei*d0,  d^,  a$,  d^.  Here  i 


17,  corresponding  to  U .25 
seconds,  roughly  one-third  of  the  mean  wave  period.  The  resulting  prediction 
may  be  compared  with  the  observed  time  history  at  the  second  point  in  Figure  23a. 


Further  Use  of  Method 

For  the  28,5-foot  ocean  data  an  eight- term  prediction  kernel  was  also  computed. 
The  improvement  in  accuracy  of  this  kernel  over  the  previous  four-term  kernel 
did  not  appear  to  be  significant  enough  to  warrant  the  additional  labor  of 
computation.  For  the  same  data  a prediction  kernel  of  a different  nature  was 
computed,  being  so  determined  that  the  squared  error  of  prediction  was  a minimum 
when  averaged,  not  over  the  entire  time  hiatory,  but  over  only  those  times  when 
the  fluctuating  surface  deviated  from  its  average  value  by  at  least 
times  its  r.m.s.  deviation  from  its  mean  value.  Since  \[Cn7a.y  X has  been 

found  in  other  studies  to  correspond  to  the  height  of  an  average  wave  measured 
from  the  level  of  its  mean  ordinate,  when  the  total  energy  in  the  spectrum  is 

K o , this  kernel  is  designed  to  be  accurate  for  higher- chan-average  waves. 
Although  in  principle  such  a kernel  would  be  more  suitable  than  the  ordinary 
least-squares  kernel  if  high  waves  are  of  especial  importance,  when  applied 
to  the  present  data  there  appeared  to  be  little  systematic  difference  in 
performance  whsn  the  waves  were  high. 


The  205-foot  ocean  wave  data  was  subjected  to  a test  of  the  assumption  of 
temporal  homogeneity,  that  is,  that  of  its  being  stationary  in  time.  In  addi- 
tion to  the  1000-second  time  histories  referred  to  above  there  was  also  measured 
at  the  seme  two  points  time  histories  of  6Q-seconds  duration  beginning  approxi- 
mately 30  minutes  after  the  end  of  the  earlier  ones.  No  estimate  was  made  of 
the  correlation  functions  for  these  later  time  histories,  but  the  prediction 
kernel  commuted  for  the  earlier  date  was  applied  to  the  60-second  data.  The 
results,  shown  in  Figure  23b  indicate  that  for  purposes  of  linear  least-square 
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prediction,  the  correlation  functions  have  not  changed  significantly  over  a 
30-minute  interval. 

The  possibility  of  using  the  recently-developed  Wave  Ordinate  Distribution 
Analyser*  for  analog  computation  of  estimates  of  the  correlation  function  was 
also  investigated  by  means  of  the  28.5-foot  ocean-wave  data.  These  estimates, 
based  on  the  distribution  of  sums  of  the  time-history  functions  with  various 
time  lags  between  them,  *ould  be  relatively  easy  to  compute.  The  results  indi- 
cated that  the  estimates  w-u-d  bo  accurate  enough  to  yield  prediction  kernels 
of  the  surface  time  Wr.jry  comparing  favorably  with  those  baaed  on  directly- 
computed  values  of  and  % ^ ' . 


Extensions  of  the  Method 

Modifications  of  the  leas't-squares  prediction  methods  described  thus  far  may  in 
certain  cases  be  expected  to  yield  somewhat  better  predictions.  For  example, 
one  may  carry  out  linear  prediction  of  the  time  history  at  one  point  from  observed 
time  histories  at  two  or  more  seacard  points,  a method  which  is  perhaps  suitable 
for  non-uniform  short-crested  wave  systems.  An  extension  of  the  correlation 
function  concept  leads  to  non-linear  prediction,  . method  which  may  be  suitable 
for  shallow-water  prediction. 

If  the  requirement  to  predict  individual  instantaneous  time-history  ordinates 
is  relaxed,  somewhat  better  prediction  may  be  possible.  Two  Alternative 
quantities  whose  prediction  suggests  itself  are  (1)  the  envelope  of  the  time- 
history  function  and  (2)  the  square  of  the  ordinate  integrated  over  a constant- 
length  interval  whose  endpoints  advance  with  time  and  whose  duration  is  long 
compared  with  a wave  period.  The  former  would  correspond  to  individual  waves 
in  the  wave  record  while  the  latter  would  be  proportional  to  a moving  average 
of  the  energy  over  a fixed  time  interval. 

Discussion  of  the  Method 

The  advantages  of  the  least-squares  prediction  methods  for  wave  transformation 
are  the  relative  shortness  of  the  kernels  required  for  accurate  prediction,, 
and  their  general  applicability.  They  may  be  used  on  data  taken  under  many 
conditions,  whether  capable  of  being  dealt  with  by  hydrodynamic  theory  or  „. 

High  accuracy,  when  obtainable,  may  be  attributed  to  the  fact  that  the  method 
provides  for  each  situation  a solution  which  is  tailor-made  to  it. 


The  chief  disadvantages  of  the  least-squares  method  are,  of  course,  related  to 
the  specialness  of  the  solution,  in  that  special  information  — namely,  values 
of  the  correlation  functions — must  be  known  before  it  can  be  used.  This 
information,  which  is  equivalent  to  a knowledge  of  certain  features  of  the 
spectra  of  the  two  time  histories  involved,  is  required  for  each  situation  in 
which  the  method  is  to  be  applied.  The  labor  of  computing  in  each  situation 
th9  correlation  functions  and  the  prediction  kernel  itself  is  an  inconvenience 
which  nay  be  lessened  by  the  utilization  of  analog  computers  for  the  nurpose 
(e,gt,  the  ordinate  distribution  analyzer  mentioned  above). 


Tn  those  situations  where  it  is  known  that  hydrodynamic  theory  can  predict 
the  wave  transformation  that  occurs  between  the  two  points,  this  knowledge,  . 
may  be  used  to  determine  any  desired  values  of  the  correlation  function  ^ 
provided  the  correlation  function  is  known  completely  for  the  seaward 

point.  In  such  a ease  a table  may  then  be  compiled,  once  and  for  all,  showing 
*Lun£,  W.W.,  "An  Electronic  Instrument  for  Statistical  Analysis  of  Ocean 
Waves",  University  of  California,  Institute  of  Engineering  Res ^ 
Technical  Report,  Series  3,  Issue  3h 3,  September  1952.  ti  “ 
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the  least-squares  prediction  kemol  to  be  used  for  any  given  "input"  spectrum, 
since  the  latter  determines  Since  this  spectrum  also  determines 

those  distances  and  times  over  which  prediction  will  be  most  accurate,  this 
information  may  likewise  be  tabled  in  such  a case. 

It  remains  to  obtain  a knowledge  of  the  input  spectrum  (or  correlation  function), 
a problem  which  may  be  approached  in  various  ways.  In  view  of  the  results  shown 
in  Figure  23b  an  estimation  of  the  correlation  function  based  on  data  taken  a 
short  time  (perhaps  up  to  hour)  before  the  prediction  was  to  be  made  would 
probably  be  sufficient.  If  the  estimation  of  the  correlation  function  were 
made  even  earlier,  the  change  reflected  in  the  values  occurring  later  could  un- 
doubtedly be  forecast  to  some  extent  from  a knowledge  of  the  past  values  of 
the  function  and  from  meteorological  conditions.  When  adequate  statistical 
data  on  ocean-wave  spectra  have  been  accumulated,  sufficient  regularity  in  their 
shape  may  be  found  to  permit  their  characterization  by  a few  parameters,  and 
the  rate  of  change  of  the  latter  with  time  may  prove  to  be  sufficiently  slow 
to  penult  the  prediction  of  spectra  a short  time  in  advance. 


Summary  Based  on  Preliminary  Work 

Examples  presented  above  show  that  linear  least-squares  prediction  of  time- 
history  ordinates  can  be  applied  to  wave-transformation  problems  to  yield  mod- 
erately accurate  Predictions  in  the  cases  considered.  When  prediction  over 
greater  times  and  distances  is  attempted,  the  above-mentioned  extensions  of 
the  method  may  be  needed. 

While  the  main  purpose  of  the  work  described  in  this  section  on  least^squares 
prediction  was  to  investigate  the  degree  of  predictability  inherent  in  phenomena 
of  wave  transformation,  certain  modifications  of  the  techniques  already  employed 
have  been  mentioned  in  order  to  indicate  that,  should  other  methods  of  prediction 
prove  inadequate,  the  least-squares  method  may  perhaps  be  made  serviceable.  Much 
of  the  practical  success  of  the  method  would  eeem  to  depend  on  the  accumulation 
of  knowledge  of  the  form  and  frequency  of  occurrence  of  wave  spectra  under 
various  conditions. 
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